Surface segregation in nanoparticles from first principles 
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FePt nanoparticles are known to exhibit reduced Llo order with decreasing particle size. The 
reduction in order reduces the magnetic anisotropy and the thermal stability of the direction of 
magnetization of the particle. The phenomenon is addressed by investigating the thermodynamic 
driving forces for surface segregation using a local (inhomogeneous) cluster expansion fitted to ab 
initio data which accurately represents interatomic interactions in both the bulk and surface regions. 
Subsequent Monte Carlo simulations reveal that first surface layer Pt segregation is compensated 
by Pt depletion in the second subsurface layer. This indicates that the core's ordered state is 
not affected by surface thermodynamics as much as previously thought. Thus, the weak ordering 
experimentally observed is likely not due to fundamental thermodynamic limitations but rather to 
kinetic effects. 
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Understanding the physics of the transition from bulk 
to nanoscale magnetic alloys is important from both a 
fundamental [ij and technological point of view. In the 
last decade, Fe-Pt nanoparticles have been intensively 
investigated in view of possible future applications as an 
ultra-high density information-storage medium and high- 
performance permanent magnets U [1, |3| • The critical 
issue for information-storage application is the presence 
of magnetic anisotropy offering sufhciently large thermal 
stability. In Fe-Pt, a high magnetic anisotropy is guar- 
anteed by an ordered LIq phase. However, recent exper- 
imental observations have shown a difficulty in obtaining 
a high degree of LIq order in FePt nanoparticles annealed 
at T <600°C with diameter less than ~ 4nm Due 
to the high surface to volume ratio of nanoparticles, sur- 
face segregation has been suggested to be one of the possi- 
ble causes of this reduced ordering 0, 0, [1| (interestingly, 
in FeC and FeMoC nanoparticles, size plays the opposite 
role inducing a disorder-order transition [10]). 

In this letter, we address the FePt surface segregation 
problem with accurate ah initio computational thermo- 
dynamics techniques. In earlier studies, the phenomenon 
was tackled using several interatomic model potentials 
0) B @] ■ Pair interactions were fitted to the experimen- 
tal phase diagram [7, 8], and to bulk first principles data 
01 . A single-layer isotropic surface potential was esti- 
mated from the surface energy difference between pure 
fee Fe and PtJIJ and by fitting to experimental segrega- 
tion profiles [Sj] . The Embedded Atom Method potential 
was used in Ref. More recently, surface segrega- 

tion in small Fe-Pt clusters was studied by direct com- 
parison of the energies of several different configurations 
obtained from first principles [ll| . Our approach is rad- 
ically different: the surface potential is obtained from ah 
initio calculations without fitting to experimental data 
and without a-priori assumptions about strength and 
isotropy. We combine the surface potential obtained us- 
ing the method just described and the bulk pair potential 
into a local (inhomogeneous) cluster expansion [l^ . 



enabling efficient Monte Carlo simulations to describe the 
nanoparticle's actual thermodynamic equilibrium segre- 
gation profile. The temperature and size dependencies 
of the Llo order identify the regimes when surface segre- 
gation is responsible for reduced equilibrium order with 
correspondingly low magnetic anisotropy. 
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FIG. 1: The unit cells of crystal structures representing 
atomic slabs separated by vacuum used to model the nanopar- 
ticle's facets. Each slab is obtained from an unrelaxed Llo 
crystal structure along (a) (001), (b) (100), and (c) (111) 
planes. Layers (L) and atoms are labeled with numbers. Fe 
and Pt atoms are represented as black and grey circles, re- 
spectively. 



The nanoparticle is modeled as an fcc-based "trun- 
cated octahedron" (Fig. [1]), which is typically observed 
for chemically synthesized FePt nanoparticles [H, with 
free boundary conditions. The initial atomic configura- 
tion is considered as ordered LIq with alternating Fe and 
Pt layers forming (001) fee crystal planes. This parti- 
cle has three inequivalent types of surfaces, (001), (100), 
and (111), which are addressed by using standard peri- 
odic slab geometries (Fig. [1]). The vacuum layer thick- 
ness was chosen to be approximately equal to three fee 
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FIG. 2: (Color online) Formation energies of different slabs 
with respect to ideal Llo slabs from both first principles 
(GGA) and using three cluster expansions "aniso2", "iso2", 
and "bulk" (see Tab. UJ. The slabs are identified by the la- 
bel "iji" , where i and j indicate the layers between which the 
atoms are exchanged (Fig. [T]). Two configurations are con- 
sidered: ij and ij for every choice of i and j. "Pt-up" /"Pt- 
down" mean that Pt atoms move to/from the surface of the 
corresponding perfect Lfo slab, respectively. 



lattice parameters. As well as perfectly ordered slabs, de- 
fected surfaces obtained by exchange(s) of atoms belong- 
ing to different atomic layers are considered. The ener- 
gies were calculated from first principles within the Gen- 
eralized Gradient Approximation using Projector Aug- 
mented Wave pseudopotentials, as implemented in the 
VASP package [li, Il5l |. The calculated slab energy dif- 
ferences are presented in Fig. [2l In agreement with pre- 
vious theoretical 0, S [HI EB] and experimental reports 
17| , the data indicates that there is an energetic gain for 
Pt atoms to segregate into all three considered surfaces. 
However, this finding is not sufficiently compelling to de- 
termine the thermodynamic equilibrium of the system. 
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FIG. 3; (Color online) The dependence of surface potential 
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LCE ^ 'bulk layer number in (a) (001), (b) (100), 

and (c) (111) slabs corresponding to initial first principles 
(GGA) data and three cluster expansions: "aniso2", "iso2", 
and "bulk" (see Tab. U). 



To parameterize the entropic contributions to surface 
segregation, a number of local cluster expansions (LCE) 
are constructed so that they reproduce atom exchange 
energies obtained from total energy calculations both in 
the bulk and in the surfaces. The local (or inhomoge- 
neous) nature of the cluster expansions (CE) manifests 
itself by the presence of a layer-dependent unary mix- 
ing potential V'-^^ (see Eqs. (4,7) in 6]). Such a layer 
dependence near the surface may be formally considered 
as an external surface potential applied to the surface 
atoms. The constructed LCEs differ from each other by 
the number of external layers affected by surface poten- 
tial and by the directional dependence of surface poten- 
tial. The constructed LCEs and their accuracy in sur- 
face regions are presented in Fig. [3] and Tab. H] The 
previously proposed bulk CE |6|] has to be modified only 
at two external surface layers (see Fig. ^ so that only 
the LCEs with one and two layers of surface potentials 
have to be considered. Comparing the prediction errors 
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TABLE I: Cluster expansions obtained for a description of 
nanoparticle configurational energies: "bulk" is the bulk CE 
obtained in ^ without accounting for surface effects; "isol" 
and "iso2" ("anisol" and "aniso2") designate the implemen- 
tation of isotropic (anisotropic) surface potential affecting one 
and two external layers, respectively. "lOOisol" and "100iso2" 
are the hybrid CEs for which the surface potential is the same 
for all (001) and (100) types of surfaces but is different from 
the surface potential for the (111) surface. The least-square- 
fitting (LSF) errors are those of atom exchange energies in 
two external layers (total and within each surface). 
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of different CEs, it is concluded that accounting for sur- 
face the potential within the second layer has a larger 
effect than considering the anisotropy in the surface po- 
tential. Hence, we conclude that the "iso2" LCE (with 
isotropic two-layer surface potential) represents the best 
compromise between the accuracy and complexity of cal- 
culations. Note that the "isol" surface potential (-0.33 
eV) is very similar to the corresponding surface poten- 
tial (-0.30 eV) obtained in Refs. [vj, Q with different 
approaches. 

The cluster expansion is then employed in a finite- 
temperature Monte Carlo simulation scheme While 
spanning the phase space, the LIq order is monitored by 
the "generalized" parameter rj: 



V = (niax{|77^|, \T]y\, \t]^\} 



-{|^.U%U^.|}>, (1) 



where (...) is the statistical average over the Monte Carlo 
steps, and three directional parameters r]i{i = x, y, z) are 
defined as the difference between the Pt atom concen- 
trations at odd and even crystal planes perpendicular 
to the i-th direction. The generalization is required for 
small nanoparticles, because, in the presence of a surface 
potential, we have observed that all rji can be compara- 
ble in wide temperature intervals. Generally, we can not 
neglect minjlry^;!, jry^l, jTy^j} in Eq. ([1]) as it was done Ref. 
0. In the case of isotropic order, \rjx\ = \'qy \ = \r]z\, we 
have ^ = 0. 

The results of our Monte Carlo simulations with the 
"iso2" cluster expansion are presented in Figs. [Uni The 
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FIG. 4: (Color online) The temperature dependence of the 
equilibrium (a) Llo order parameter rj and (b) Pt concentra- 
tions in two external layers and internal core for truncated 
octahedron ("TO") shape nanoparticles with different sizes 
at near equiatomic composition. In graph (a), we include the 
data obtained within the analytical ring approximation [3l 
for bulk ("bulk") as well as by Monte Carlo simulation for 
the parallelepiped ( "ppd" ) sample containing A*' — 60'^ atoms. 
The solid and dashed lines represent the total and core (ex- 
cluding two external layers) order parameters, respectively. 



comparison of the present and previous [6(] results con- 
firms that the presence of a surface potential reduces the 
total order in agreement with Refs. [3, [§]■ The reduction 
is larger for smaller particles having a bigger fraction of 
surface versus volume. In addition, there is a negligible 
effect of the surface potential on the core order parameter 
for particles larger than '^3.4nm, followed by a still con- 
siderable total order parameter (^^0.4-0.5 of maximum 
value at 600°C). Only for particles smaller than ~3.4nm, 
is there a strong reduction in order. 

Figure[4l^b) shows that the Pt-segregation into the first 
external layer is compensated by Pt-depletion in the sec- 
ond layer so that the composition of the core remains 
close to ideal. This is a direct consequence of the values 
of the surface potential "iso2" for the two external layers. 
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FIG. 5: (Color online). Equilibrium LIq order parameters rj 
versus particle size at T =600°C: total, core (excluding two 
external layers) and obtained in Ref. @] by neglecting surface 
potential. 



similar in magnitude but with opposite sign (Fig. [3] and 
Tab. IJ). Consequently, as the core remains close to sto- 
ichiometry, its ordered state is affected by surface ther- 
modynamics less than that previously concluded 0, Q . 

Our model considers only configurational entropy and 
neglects other effects. For instance, the effects of vibra- 
tional [3| or magnetic entropies are not included. 
However, the good agreement between the measured and 
calculated bulk transition temperatures (Fig. |4]) sug- 
gests that these corrections are small in the important 
range T < 600°C. Furthermore, we have also consid- 
ered nanoparticles embedded in vacuum rather than in 
the polymeric medium typically used. Nevertheless, our 
surface potential is similar to previous estimates based 
on experimental data Q and, more importantly, our 
most crucial observation pertains to the second subsur- 
face atomic layer, which is clearly less sensitive to the 
surrounding media. Therefore, our findings reopen the 
possibility that the weak ordering observed experimen- 
tally is not due to fundamental thermodynamic limita- 
tions but probably to kinetic effects 2]J that may be 
easier to control [6|. 
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